%script que simula el avance de temperatura de la habitacion
%sin ningun artefato de control
%aparece la temperatura de referencia solo como algo extra, pero no incide
%en la temperatura de la habitacion

Ti = 20; %temperatura interna comienza en 20�
Te = 20; %temperatura externa comienza en 20�
Tr = 20; %temperatura referencia comienza en 20�
T = 10; %periodo de muestreo [s]
tFinal = 1*60*60; %duracion de ¿seis? horas, medido en [s]
prob_puerta = 1/(36*T); %probabilidad de apertura de la puerta
n = 1; %indice
time = [0]; %vector de tiempos
for t=T:T:tFinal
    n = n+1;
    Te = [Te, externalTemperature(t)];
    Tr = [Tr, referenceTemperature(t)];
    a = randi(1/prob_puerta); % numero entero entre [1,1/prob_puerta
    if(mod(a,1/prob_puerta)==0) %la puerta esta abierta
        newT = 0.169*Ti(n-1) + 0.831*Te(n);
        disp('Se abrió la puerta');
    else
        newT = 0.912*Ti(n-1) + 0.088*Te(n);
    end
    Ti = [Ti,newT];
    time = [time,t];
end
figure(1);
plot(time,Tr);
figure(2);
plotyy(time,Te,time,Ti)
legend('Temperatura Externa','Temperatura Interna');